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Abstract 

Consider time-harmonic electromagnetic wave scattering from a biperiodic dielectric 
structure mounted on a perfectly conducting plate in three dimensions. Given that unique- 
ness of solution holds, existence of solution follows from a well-known Fredholm framework 
for the variational formulation of the problem in a suitable Sobolev space. In this paper, we 
derive a Rellich identity for a solution to this variational problem under suitable smooth- 
ness conditions on the material parameter. Under additional non-trapping assumptions on 
the material parameter, this identity allows us to establish uniqueness of solution for all 
positive wave numbers. 

1 Introduction 

Scattering of electromagnetic waves from periodic structures is not only an interesting mathe- 
matical topic in its own right but also of great interest in applications, e.g., for the construction 
and optimization of optical filters, lenses, and beam-splitters in optics. An overview about this 
and further topics in applied mathematics related to wave propagation in periodic structures 
can be found in, e.g., [5]- In this paper we consider scattering of time-harmonic electromag- 
netic waves from a dielectric biperiodic structure mounted on a perfectly conducting plate in 
three dimensions. By biperiodic, we mean that the structure is periodic in the, say, x\- and 
^-direction, while it is bounded in the x% direction. In contrast to scattering from bounded 
structures, uniqueness of solution for this scattering problem does in general not hold for all 
positive wave numbers. Instead, non-trivial solutions to the homogeneous problem might exist 
for a discrete set of exceptional wave numbers, and these solutions turn out to be exponentially 
localized surface waves. Our study in the present paper focuses on conditions guaranteeing 
the well-posedness of the full three-dimensional electromagnetic scattering problem mentioned 
above. We establish non-trapping and smoothness conditions on the (non-absorbing) dielec- 
tric such that uniqueness of solution holds for all positive wave numbers. This means that 
materials satisfying the latter conditions cannot guide surface waves. 
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Figure 1: Sketch of the biperiodic structure under consideration. 

Mathematical formulation for the well-posedness of electromagnetic scattering problem for 
periodic structures has been an active area of research in the last years. For the scalar case, 
the authors in [2j[T7] studied uniqueness of solution for all wave numbers (or, equivalently, all 
frequencies), under geometrical conditions on the scatterer, for impenetrable structures with 
Dirichlet and Neumann conditions. Similar results are obtained in the paper [7] for more 
complicated periodic structures which are constituted of conducting and dielectric materials. 
The latter paper further gave examples of structures for which non-uniqueness of solution 
occurs at the so-called singular wave numbers. These wave numbers were shown to be related 
to guided waves (surface waves) that are exponentially localized along the structure. 

For the case of Maxwell's equations, the authors of [10] studied well-posedness of the 
scattering problem from a medium consisting of two homogeneous materials separated by a 
smooth biperiodic surface using an integral equation approach. In [Hl6llll] the authors studied 
existence and uniqueness of solution for the scattering problem from penetrable biperiodic 
structures using a variational approach for the magnetic field. Nevertheless, unlike the scalar 
cases, the uniqueness results in cited cases of Maxwell's equations were proven for all but 
possibly a discrete set of wave numbers. Furthermore, all the cited papers above considered 
the non-magnetic case, i.e, the coefficient magnetic permeability is assumed to be the same 
constant outside and inside the structure. The case of variable magnetic permeability were 
investigated in the paper [1] for Maxwell's equations where the biperiodic structure consists 
of conducting and dielectric materials. That paper studied a variational approach, formulated 
in terms of the electric field, and showed that the obtained saddle point problem satisfies 
the Fredholm alternative, and again uniqueness of solution was proven for all but possibly 
a discrete set of wave numbers. More recently, the paper [16] analyzed the well-posedness 
of the scattering problem for penetrable anisotropic biperiodic structures with a restriction 
on the non-magnetic case again. The latter paper also proved that the scattering problem is 
uniquely solvable for all wave numbers if the structure contains absorbing materials, and if the 
dielectric tensor is piecewise analytic. Hence, to the best of our knowledge, uniqueness results 
for all wave numbers for the vectorial scattering problem still remains open if the biperiodic 
materials is non-absorbing. 

The aim of the present work is to prove that the electromagnetic scattering problem for non- 
absorbing biperiodic dielectric structures mounted on a perfectly conducting plate is uniquely 
solvable for all positive wave numbers if the material parameter satisfies non-trapping and 
smoothness conditions. We formulate the Maxwell's equations variationally in terms of the 



magnetic field in a suitable Sobolev space. We further restrict ourselves to the case of non- 
magnetic and isotropic materials. The variational problem is well-known to fit into a Fredholm 
framework, see, e.g., [HHTJITB]. (These papers deal with periodic scattering in the full space, 
but can be adapted to the half-space setting that we consider here.) As mentioned in the 
paper [7] on the corresponding scalar scattering problems, non- uniqueness phenomena indeed 
arise at certain singular wave numbers if the non-absorbing material parameter satisfies suitable 
trapping conditions. In this paper we show a converse result for the full three-dimensional 
periodic Maxwell equations: uniqueness of solution holds for all positive wave numbers if 
the material parameter is non-absorbing and satisfies suitable non-trapping and smoothness 
conditions. To prove the uniqueness result we derive a so-called Rellich identity for a solution 
to the homogeneous variational problem. The solution estimates resulting from this integral 
identity allow us to show that the homogeneous variational problem has only the trivial solution 
for all positive wave numbers. 

Our analysis extends the approach in [12] that was motivated by an existence and unique- 
ness proof for solutions to rough surface scattering problems via Rellich identities in [8]. For 
scalar periodic problems, a related technique has been used in [7J. The paper [12] studied elec- 
tromagnetic scattering from rough, unbounded penetrable layers. Such scattering problems 
are considered to be more complicated than those for periodic structures since the problem to 
find the scattered field cannot reduced, e.g., to a bounded domain. The applications of rough 
scattering problems include for instance outdoor noise propagation, oceanography or even op- 
tical technologies when the dielectric lacks periodicity. The authors in [T2] formulated the 
latter scattering problem in terms of the electric field. We will instead choose a formulation in 
terms of the magnetic field, which somewhat changes the role of the dielectric material param- 
eter in the integral identities since the material is non-magnetic. The paper |12] establishes 
existence and uniqueness of solution under non-trapping and smoothness conditions on the 
material parameter. While a priori estimates resulting from the Rellich identity allowed the 
authors in [12] to deduce uniqueness of solution, existence of solution has been obtained using 
a limiting absorption argument. The approach studied in the present paper is, from the tech- 
nical point of view, somewhat similar to the one introduced in [12]. However, the analysis of 
the biperiodic case is definitely simpler since uniqueness of solution directly implies existence. 
Therefore, one only needs to investigate the Rellich identity and estimates for solutions to the 
homogeneous problem. It turns out also that this procedure produces weaker assumptions on 
the material parameter than those found in [12]. More precisely, uniqueness and existence of 
solution for all wave numbers are obtained under the following (non-trapping and smoothness) 
assumptions on the biperiodic relative material parameter e T : R 3 , := {x G M 3 , x% > 0} — > M.. 
First, we assume that e~ l € L°°(1R^_) equals one in {x$ > h} for some h > and possesses 
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essentially bounded and measurable first weak derivatives. Second, we require that 
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where V^e^ 1 := (de~ l /dxi,ds~ 1 /dx2,0) T ■ Under these conditions, the existence of surface 
waves is automatically ruled out. While conditions (a) and (c) are similar to conditions (a) 
and (d) in [121 Eq. (7.2)], condition (b) is weaker and clearly simpler than the corresponding 
conditions (b) and (c) in [T2J Eq. (7.2)]. 

The half-space setting that we consider in this paper is somewhat special, and it seems 
worth to mention that the Rellich identity itself generalizes to a corresponding periodic scat- 
tering problem in full space. The resulting estimate for a solution H to the scattering problem 
has a similar structure to the estimate in Lemma 15.41 However, in the half-space setting, the 
term 2Re J^{de~ l / dx^) (dH^/dx^H^ dx can be treated without integration by parts using a 
Poincare lemma. In contrast, in the full-space setting the only obvious way of treating this 
term is to integrate by parts. Since we seek for solution estimates, this introduces the condition 
that £3 1— > e~ 1 (x±, X2, x$) needs to be concave to conclude. Since this is a somewhat unnatural 
condition, we do not present this result in more detail. 

One can further generalize the results presented here to certain anisotropic structures. 
However, already for the simpler case of isotropic coefficients the derivation of the Rellich 
identity is a technical matter. Again, we have opted to try to keep the presentation simple 
instead of treating the most general setting that could be considered. 

The paper is organized as follows: In Section 2 we present setting of the problem. Section 

3 is dedicated to a variational formulation and to the Fredholm property of the latter. Section 

4 contains a couple of technical lemmas. We derive the integral inequalities resulting from the 
Rellich identity in Section 5. Finally, the uniqueness of the variational problem for all wave 
numbers is proven in Section 6. 

Notation: We denote by H s (R d ) 3 , d = 2,3, the usual L 2 -based Sobolev space of vector- 
valued functions in R d . Moreover, Hf oc (R 3 ) 3 = {v E H S (B) 3 for all balls B C M 3 }, and 
Wi.°°(R3) = { v e L°°(M. 3 ) : Vv 6 L°°(IR 3 ) 3 }. 



2 Problem Setting 

We consider scattering of time-harmonic electromagnetic waves from a biperiodic structure 
which models a dielectric layer mounted on a perfectly conducting plate. The electric field E 
and the magnetic field H are governed by the time-harmonic Maxwell equations at frequency 
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tu > in R^ = {(xi,x 2 ,x 3 ) G R 3 : x 3 > 0}, 

curlH + iueE = in R^, (1) 
cwlE -iufiH = Q inR%, (2) 
e 3 xE = on {x 3 = 0}, (3) 

where e 3 = (0, 0, 1) T . The electric permittivity e is a bounded measurable function that is 2tt- 
periodic in x\ and x 2 . Further, we assume that e equals Eq > outside the biperiodic structure, 
that is, for x 3 > h where h > is chosen larger than sup{x3 : (x±, x 2 , x 3 ) T G supp(e — £o)}. 
The magnetic permeability = fiQ is assumed to be a positive constant and the conductivity is 
assumed to vanish. As usual, the problem (P)-© has to be completed by a radiation condition 
that we set up using Fourier series. 

The biperiodic structure is illuminated by an electromagnetic plane wave with wave vector 
d = (di,d 2 ,d 3 ) 6 R 3 , d 3 < 0, such that d • d = uj 2 eQfiQ. The polarizations p, q 6 R 3 of the 
incident wave satisfy p ■ d = and q = l/(a;eo)(p x d). With these definitions, the incident 
plane waves E % and H % are given by 

E i := qe ld - x , £P := pe id ' x x G R 3 . 

In the following we will exploit that one can explicitly compute the corresponding reflected field 
at {x 3 = 0}. To this end, we introduce the notation a = (01,02, — as) T for a = (01, a 2 , a 3 ) T G 
R 3 . The reflected waves at the plane {X3 = 0} are 

E r (x) := -qe llx , H r {x) := pe ilx , x G R+ , 

since divE r = 0, divH r = 0, and e 3 x (E i + E r ) = 0, e 3 • (H i + H r ) = on {2:3 = 0}. From 
now on, we denote the sum of the incident and reflected plane waves by 

E ir := E l + E r and H ir := H l + H r . 



Set 

q = (ai,a 2 ,a 3 ) T := (di,d 2 ,0) T 

and define E l £ and by 
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such that E^ and H£ are 27r-periodic in x\ and x 2 . If we apply the same phase shift to 
solutions E and H of the Maxwell equations (HI)-©, 

E a = e- iax E(x), H a = e- iax H(x), 

and if we denote 



V Q / = V/ + iaf, curl Q F = curl F + ia x F, div Q F = divF + ia ■ F 
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for scalar functions / and vector fields F, then E a and H a satisfy 

curl Q H a + iueE a = in W\_, (4) 



curl Q E a — ioofiQH a = in (5) 
e 3 x E a = on {x 3 = 0}. (6) 

Note that we still have div a curl Q = and curl a V a = 0. Let us denote the relative material 
parameter by 

£ 

£ r := — . 

Obviously, e r equals one outside the biperiodic dielectric structure. Recall that the magnetic 
permeability jjlq is constant which motivates us to work with the divergence-free magnetic field, 
that is, drv a H a = 0. 

Note that @) plugged in into (0 implies that e 3 x (e" 1 curl a H a ) = on {X3 = 0} and 
that the condition e 3 ■ H a = on {^3 = 0} can be derived by plugging © into ©. Hence, 
introducing the wave number k = oj(eofio) 1 ^ 2 , and eliminating the electric field E a from (HD-dHJ), 
we find that 

curl Q (e^ 1 curl Q H a ) - k 2 H a = in R s + , (7) 
e 3 x (s~ l curia H a ) = on {x 3 = 0}, (8) 
e 3 • H a = on {x 3 = 0}. (9) 

We wish to reformulate the last three equations in terms of the scattered field H a , defined by 
H a := H a — Since, by construction, curl Q curl Q H% — k 2 H£ = in M+, H£ ■ e 3 = and 
e 3 x (e" 1 curl a H^) = on {2:3 = 0}, a simple computation shows that 

curLj (sjT 1 curl a H a ) - k 2 H s a = - curl a ((e^ 1 - 1) curl a H%) in R\, 

e 3 x (a" 1 curl Q H s a ) = on {x 3 = 0}, (10) 
e 3 • H° a = on {x 3 = 0}. 



Due to the biperiodicity of the right-hand side and of e r , we seek for a biperiodic solution H a , 
and reduce the problem to the domain (0, 2tt) 2 x (0,oo). We complement this boundary value 
problem by a radiation condition, see also in [SHU], that we set up using Fourier series. The 
scattered field H a is 2-7r-periodic in x\ and xi and can hence be expanded as 

Hi(x) = H n (x 3 )e in - X , x = ( Xl ,x 2 ,x 3 ) T G R% , A = I 2 x {0}, (11) 

n€A 

where the Fourier coefficients H n (x 3 ) are defined by 

Hn{x 3 ) = ^J J H^(x 1 ,x 2 ,x 3 )e- in - x dx 1 dx 2 , n G A. (12) 

Define 



I Jk 2 - In + a\ 2 , k 2 > \n + a\ 2 , 

Pn ■= { V r , ' ' ~ ' n G A. 

I !Y|n + ap-K , Ar < |n + a| , 
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Figure 2: Geometric setting for electromagnetic scattering problem from a biperiodic dielectric 
structure mounted on a perfectly conducting plate (in two dimensions, for simplicity). 

Since e" 1 equals one for x 3 > h it holds that div a H^ vanishes for £3 > h, and equation (llpp 
becomes (A a + k 2 )H^ = in {2:3 > h}, where A Q = A + 2ia ■ V — \a\ 2 . Using separation of 
variables, and choosing the upward propagating solution, we set up a radiation condition in 
form of a Rayleigh expansion condition, prescribing that iJ® can be written as 

H s a (x) = ^ H ne iM*a-h)+in-x for { X3 > ^ where ^ . = H n (h), (13) 

and that the series converges uniformly in compact subsets of {#3 > h}. 

The scattering problem to find a scattered field that satisfies the boundary value 
problem (|10p and the expansion (|13p is in the following section reformulated variationally in 
a suitable Sobolev space. 

3 Variational Formulation 

We solve the scattering problem presented in the last section variationally, and briefly recall 
in this section a variational formulation of the problem in a suitable Sobolev space. Our 
framework is an adaption of the results from [16] to our half-space setting. In contrast to the 
variational formulation in //(curl) in pQ, the papers [^l6"llll|[T6] set up a variational formulation 
in H 1 for the magnetic field. Indeed, since the latter is divergence-free, any solution that is 
locally //(curl) indeed belongs locally to H . For our purposes, the H 1 formulation has the 
additional advantage that it is well-defined at Rayleigh- Wood frequencies, as it was noted 
in [JB]. We define a bounded domain 

Q = (0,2-7r) 2 x (0, h) for h > supja^ : (xi,X2,X3) T £ supp(e r — 1)}, 

with boundaries Tq := (0, 2ir) 2 x {0} and Yh := (0, 2ir) 2 x {/i}, and Sobolev spaces 

Hp(Q) 3 := {F E H\nf : F = F\q for some 2^-biperiodic F G H( oc (R 3 ) 3 }, t G N, 
Hl T (nf := {F = (F l ,F 2 ,F 3 ) T G //i(O) 3 : F 3 = on r }, 

equipped with the usual integral norm, e.g., 

11-^11^1(0)3 = II^IlL 2 ^)^ + 11^-^1112(0)3. 
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The space i?p T (Q) 3 of periodic vector fields that are tangential on To is well-defined due to 
the standard trace theorem in H . We also define periodic Sobolev spaces of functions with 
d = 1, 2, 3 components on IV for s£l, 

H s v (T h ) d := {F G H s (T h ) d : F = F\ Fh for some 2vr-biperiodic F G Hf oc ({x 3 = h}) d }. 

A periodic vector field F G H s (Th) d can be developed in a Fourier series, F(x) = 
EneA-^expfm ■ x), and \\F\\ H s {rh) d = (£ neA (l + n 2 ) s |^„| 2 ) 1/2 defines a norm on H°(T h ) d . 
We define a non-local boundary operator T a (the exterior Dirichlet-Neumann operator) by 

(T a f)(x) = iPnf n e in - x , for / = / n exp(m • x) G ^ /2 (T h ). 

1/2 — 1/2 ■ 

It is a classical result that T a is bounded from H p (T/j) into -ff p (T^), see, e.g., [3]. Using 
T Q , we define a vector of (pseudo-) differential operators R a := (d a /dx\, d a /dx2, T a ). For a 
vector field F G Fp /2 (r h ) 3 , 

R a xF = (d a /dx 1 ,d a /dx 2 ,T a ) x F, R a ■ F = {d a /d Xl , d a /dx 2 , T a ) • F. 

Since all components of R a are bounded operators from Hp^iT^) into H p 1 ^ 2 (Th), the operator 
F H> x F is bounded from H^ 2 (Th) 3 into i/ p " 1//2 (r/ l ) 3 , and F ^ R a ■ F is bounded 
from Hp^ 2 (Th) 3 into H p l ^ 2 (T} l ). If a biperiodic function H G //^(M^ ) satisfies the Rayleigh 
expansion condition, then T a FLj, = dH-^/dx^, on T^. This implies that e 3 x (curl Q H) = 
e 3 x (R a x F) on T h (see, e.g., [JB]). 

Assume that H a is a distributional periodic solution to the boundary value problem ([TO]) 
such that H a , curl a ii^, and div a H a are locally square-integrable, such that the radiation 
condition (fTB"]) is satisfied, and such that v ■ (H a + H%) and v x (e" 1 zm\{R a + fl£ r )) are 
continuous over interfaces with normal vector v where e r jumps. As noted in [16] . this implies 
that, following the above notation, H a G i?p T (f2). Then the Stokes formula [HQS] implies 
that 

/ (e~ l curl Q H s a ■ curl Q F - k 2 H s a ■ F) dx 
Jn 

- / e 3 x (ejT 1 cml a H s a ) ■ Fdx + / e 3 x (R a x H s a ) ■ F ds 

= [ (1 - e~ l ) curl Q JZ*" • curl a Fdx - / (e 3 x (1 - ^T 1 ) curl a • Fdx 
Jn Jr 

for all test functions F G -ffp T (f2) 3 . Since we assumed that 

= e 3 x (e" 1 curl Q H a ) = e 3 x (e^ 1 curl Q (iI* + H%)) on T , 
the above identity simplifies to 

f (e~ l curl Q H a ■ cm\ a F - k 2 H s a ■ F) dx + f e 3 x (R a x H s a ) ■ F ds 
Jn Jv h 

= [ (1 - e~ l ) curl a fl£ • cur l a Fdx - [ (e 3 x curl a H%) ■ Fdx . 
Jn Jrn 
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By construction, e 3 x curl a vanishes on To, that is, we can neglect the last term in the last 
equation. The divergence constraint div a H^, = that follows from ()10p shows that 



B(H s a , F) := [ (e' 1 curl Q H* ■ curl a F - k 2 H s a -F)dx + p [ (div a £T£)(div a F) dx 
Jn Jn 

+ I e 3 x(R a x H s a ) • Fds - [ (R a • H s a ){e 3 ■ F) ds 

= [ (1 - e' 1 ) curl a H% ■ curia F dx , (14) 
in 

where p is some complex constant with Re (p) > c > and Im (p) < 0. 

We next prove that the bounded sesquilinear form B : T (f2) 3 x x(^) 3 ~~ ^ ^ satisfies 
a Garding inequality (this goes back to pQ), i.e. there exist strictly positive constants c\ and 
C2 such that 



Re (B(H, H)) > ci / |V Q #| 2 dx-c 2 / |F| 2 dx. 

in in 



(15) 



for all F G i^ T (ft) 3 . 

Theorem 3.1. Assume that e" 1 G L°°(S1) is positive and bounded away from zero. Set 
Re p = info e" 1 > and choose Im p < 0. T/ien i3 satisfies (|15p . 

Proof. As in |16|, proof of Theorem 1] one shows that 

Re (£(#,#)) > Re(p) / (| curl a ff| 2 + |div a i/| 2 ) dx - k 2 [ \H\ 2 dx 

Jn Jn 

r r / act tt act tt \ 

-Re / T a H-Hds-2Re (H~ 3 — -± + H~ 3 — -2 ) ds . 

Jr h Jr h V dx i dx 2 J 

The following identity follows from integrations by parts, the periodicity, and the vanishing 
normal component of H on To, 

/ (| curl Q F| 2 + |div Q #| 2 )dx = I \V a H\ 2 dx +2Re j \W 3 ^^ +1T 3 ^^\ ds . 

In consequence, 

Re (#(#,#)) > Re(p) / |V Q #| 2 dx-£; 2 / |i^| 2 dx 

in in 



Re / T.ff • H is - 2(1 - Re (,)) Re / (^i + ^) ff 3 d. . 



Precisely as in |16] one shows now by a Fourier series argument that 



Re/ T a H-Hds -2(1 -Re(p))Re / ( + ^-^-) H 3 ds > Re [ K(H)-Hds 

Jr h Jr h V Ssi / ir h 



> -C / Itfl 2 dx 
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for a finite-dimensional operator K on ifj^/ 2 (T/j) 3 . Note that the last inequality follows from 
| J T K(H) • H ds | < C \H\ 2 dx due to the finite-dimensional range of K and the fact that 
on finite-dimensional spaces all norms are equivalent. The last inequality implies a Garding 
inequality for B. □ 

For simplicity we write from now on H for the searched-for scattered field in (|14p and 
replace the source function curliJ^" by a G £ H^(fl) 3 . The last theorem implies the following 
corollary. 

Corollary 3.2. The variational problem to find H £ H^^Q) 3 such that 

B(H, F)= [ (1 - eJ x )G ■ curia F dx for all F £ Hi T (fi) 3 (16) 

satisfies the Fredholm alternative, i.e., uniqueness of solution implies existence of solution. 

Note that this formulation corresponds to the usual variational formulation of the Maxwell 
equations with perfectly conducting magnetic boundary conditions in smooth bounded do- 
mains, see, e.g., [21 Section 4.5(b)]. For special material parameters e~ m 

W p 1,cc (fi) := {/ £ L°°{Sl) : f = /|n for some 2^-biperiodic / £ W 1,0 °(R 3 )} 

we will in the sequel of the paper establish a uniqueness result via a Rellich identity. The next 
lemma will be useful when proving this identity. 

Lemma 3.3. Assume that e" 1 £ Wp'°°(f2) is positive and bounded away from zero, and that 
G £ Hp(Q) 3 . Then a solution H £ i?p T (f2) 3 to problem (|16|) satisfies 

curl Q (e7 1 curl a #) - fc 2 # = curl Q ((l - e^G) in L 2 (0) 3 , (17) 
div Q # = inL 2 (n), (18) 
e 3 x (e^ 1 curl Q H) = e 3 x ((1 - e^G) i n H^ 2 (r ) 3 , (19) 
e 3 • = in ^ /2 (r ). (20) 

Moreover, 

e 3 x R a x H = e 3 x curl a in H~ 1/2 (r h ) 3 and R a ■ H = in H~ l / 2 {T h ), (21) 
and dH/dx 3 = T a (H) holds in Hp 1/2 (T h ). 

Proof. The proof that div a H = is analogous to the proof of |16( Theorem 2]. In consequence, 
using a test function F £ Cq°(£1) 3 in the variational problem (|16p shows that the solution H 
satisfies the differential equation (|17j) in the distributional sense. Since H £ H^ T (Q) 3 , (|17j) 

holds in the L 2 -sense if the right-hand side belongs to L 2 (0) 3 , which holds if e~ l £ Wp'°°(Q) 
and G £ H^nf. 
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Multiplying (|17|) by F 6 H'z T (S1) 3 , using the Stokes formula, and subtracting the resulting 
expression from the variational formulation (|16p . we find that 



/ e 3 x(R a x H)-Fds - (R a ■ H) (e 3 • F) ds - / e 3 x curl Q if • F ds 
Jr h Jr h Jr h 

+ / e 3 x (e" 1 curl a H) ■ Fds - e 3 X ((1 - e^)G) ■ F ds = 0. 

If we choose F such that F\^ = 0, then we see that e 3 x (e" 1 curia H — (1 — e~ 1 )G) = 

in H~ 1/2 (r ). If e 3 • F| rft = 0, it follows that e 3 x (R a x H) = e 3 x curl a F in H~ 1/2 (T h ) 3 . 

Hence, f? Q • H = in i/ p (I\)- These identities imply that dH/dx 3 = T a (H) in if p (r^) 
due to [TBJ Lemma 1]. □ 

Remark 3.4. Instead of the above variational formulation in ffp T (0) 3 , one can also consider 
formulations in 7/ p (curl a , fi) 3 , the natural energy space for the second-order Maxwell equa- 
tions ()10p . see, e.g., |7j/. In if p (curl Q , Q) 3 there is no bounded trace operator for the normal 
component of the field, and in consequence, the formulation (|16p needs to be adapted. Usually, 
one replaces F \— > e 3 x (R a x F) x e 3 by Q(e 3 x H), where Q is a bounded operator between 
the natural trace spaces H ' ^(r^) and H l^i^h), defined by 

(QF)(x) = -J2 ^-i k2p T,n ~ [(n + a) ■ F n ](n + a)}e inx , for F(x) = £ F n e in ' x , (22) 

nGA %Pn neA 

see, e.g., JI]/. Obviously this definition only makes sense if all f3 n are non-zero. If this is the 
case, then the variational formulation (|16j) is equivalent to the formulation in H p (cuvl a , J7) 3 
obtained using Q. Under the assumption that f3 n ^ 0, all subsequent results could also be 
obtained via the formulation in i/ p (curl a , O) 3 . 

4 Integral Identities 

This section is concerned with technical lemmas that will be used to derive the Rellich identity 
and solution bounds subsequently. Roughly speaking, for deriving the Rellich identity, we 
will multiply the Maxwell equations (|17p by x 3 dH/dx 3 and integrate by parts. Therefore, 
it is the aim of the technical lemmas in this section to analyze the term Re J^x 3 dH/dx 3 ■ 

curl Q (e^ 1 curia H) dx for a solution H € T (f2) 3 to the problem (|16p. Note that the first two 
lemmas need the function H to be in H^(Q) 5 . These lemmas for the magnetic field formulation 
actually correspond to the ones for the electric field formulation in [12^ Section 3]. 

We need to introduce some notation. For a vector field F = (iq, F2, F 3 ) T we denote by 
Ft = (-Fi,i<2,0) T its transverse part. Recall that d a f/dxj = df/dxj + ictjf for a scalar 
function / and j = 1,2,3. Further, we introduce 

/ df df \ T (d a f d a f \ T — > fd a f d a f \ T 

\ax\ 0x2 J \ox\ 0x2 ) \ox2 ox\ ) 
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and, for a vector field F = (Fi, F2, F$) T 7 



d a Fi d a F2 d a F2 d a F\ 

div Qi rF := — 1 — - — and curler F := 



dx\ 8x2 



It is straightforwards to show that diVc^rcurLj,^ = as well as curl^T V Qj r = 0. Moreover, a 
tedious computation shows that 



curl Q F = (curl QjT F T )e 3 + curl a>T -F3 

and further 

I curl Q F\ 2 = I curl Q t Ft\ 2 + |curl Q tF^\ 2 + 



djF x e 3 ) 
dx 3 



dF T 



dx? 



2Re [ V a rF 3 



dF T 
dx 3 



(23) 



Lemma 4.1. Assume that e r 1 G Wp'°°(£l) is positive and bounded away from zero and that 



H G H^n) 3 . Then 



dH 



2Rc / x 3 — — • curl Q (s T 1 curl Q H) dx 
0x3 



d{x 3 e T x ) , p H 2 dx + h j I cmla H , 2 dg 



dx? 



L 



+ 2Re I e~ l ( e 3 x • curl Q H dx + 2/iRe f • (e 3 x curl a H) ds . (24) 

J si V 9xsJ Jr h dx 3 

Proof. Denote by v the outward unit normal to f2. Using integration by parts and noting that 
v = es on Th, and that the boundary term on Tq vanishes since X3 = on Tq, we find that 



r dH 

2Re / X3— — • curl Q (e7 1 curl a H) dx 
Jn dx 3 

= 2Re / e~ l curl Q I X3 — — I • curl a H dx + 2Re / 
Jn V 9x3 J J a 



dH 



X3— — • (y x e r 1 curl Q H) ds 

an 0x3 



_i d\cuil a H\ 2 

£ r X 3 ~ 

OX3 

+ 2/iRe 

d{x3E~ l 



dx + 2Re / e r 1 ( e% x ] • curl Q H dx 
Jn V 9x 3 J 



dH T 

9X3 



(es x curb, H) ds 



9X3 



curl a H | 2 dx + 2Re / e r 1 ( e3 x - — ] • curl Q H dx 



dH 



dx-~ 



h / I curia flrds + 2hRe 



dH T 

dx3 



(es x curl a H) ds . 



□ 



The next lemma continues the analysis of the term Re JqS t 1 {e 3 x dH/dx%) • curl a ffdx 
in the right hand side of 
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Lemma 4.2. Assume that e r 1 £ Wp'°°(f2) is positive and bounded away from zero. Then for 
all H S i?p(0) 3 the following identity holds, 



2Re / e r 1 ( e 3 x ] • cmi Q H dx = 2 / e r 1 
A} V 9x 3 / J n 

d(e- l Th) 



dH 



-2Re 



<9x ? 



<9x 3 

div a H dx — 2Re 



dx + 2Re / Ve" 1 • — #3 dx 
dx 3 



n 

dH 3 



div a Hj H 3 ds 



-2Re / e^Hsdiv^THT ds . (25) 



Proof. First, we have 



2Re 



Jn 



e 3 x - — I • curl a H dx = 2 / e r 1 
9x 3 / 



<9x ? 



dx 



-2Re / e. 



_! dH T — 



3X3 



• V T # 3 dx + 2Re / e 



-i^r ,.. TF - 

3X3 



• iai^3 dx . 



(26) 



Second, we compute that 



2Re / e 



H 3 dx 



2Re / e r _1 div T 



V 9x 3 J 



H 3 dx + 2Re / V T eZ 



-1 dtf T _ 



<9x? 



F 3 dx 



-2Re 



+ 2Re / Vte, 



<9x 



-1&H3 



H^divrHT dx — 2Re / e r — — divT-f^T dx 



9xs 



-1 dH T —^ 



dx? 



H 3 dx+2Re / H 3 dw T H T ds - 2Re / e^^adivriards 

r 



Now, using the identity div^-Hr = —dH 3 /dx 3 + div a ff — ia ■ H, we obtain that 



2Re / e^ 1 ^- 1 • V T H 3 dx = 2Re 



S2 



2Re 



9x a 



—^H 3 (ia-H)dx +2Re 

3X3 



de- 1 0H 3 

~a H 3~x — dx 

dx 3 dx 3 



dx? 



xdH 3 



r —H 3 div a Hdx+2Re / e~ L — -(ia ■ H) dx + 2Re / eZ 

ox? 



dH 3 



dx? 



dx 



dH 

e^-^divaHdx + 2Re 



1 . ^^H 3 dx + 2Re H 3 dw T H T dx 

dx 3 



2Rc / £Z~ L H 3 div T H T ds 
Jt 
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Applying Green formula to the term 2Re f^(de r 1 /dx^H^iia ■ H) dx , we have 



2Re 
- 2Re 



0x3 



-2Re 



de^ 1 

— —H 3 div a Hdx -2Re 

ox* 



2Rc 



n 
<9x 3 



e r — — • iaHz dx + 2Re 
dx 3 



dm 



8x3 



dx 



div a H dx + 2Re 

0x3 

div a H)J{^ds - 2Re 



r 



Ve 7 l ~lhdx 
0x3 

e7 Hsdiva^HT ds 



Now the claim follows from substituting this identity into equation (|26p . □ 

In the following final lemma of this section we will reformulate the term Re J n x^dH/dx^ ■ 

curl Cf (e^ 1 curio, H) dx for a solution H G .£fp T (f2) 3 to the problem (|16p using the last two 
lemmas. 

Lemma 4.3. Assume that e~ l G Wp ,oc (Jl) is positive and bounded away from zero. Then any 
solution H G //p T (il) 3 to i/ie problem (|16p satisfies 



2Re 



X3— — • curl a (e r 1 curio, if) dx 
dx 3 



d(x 3 e r ^ 
8x3 



cux\ a H\ 2 dx + h I I curio H | 2 ds 



+2 



<9# 



<9xs 



dx+2Re / Ve r _1 • — #3 dx - 2Re / i^— - ds 

Jn dx 3 J Vh 9x 3 



+2/iRe 



dH T 

8x3 



(e3 x curio H ) ds . 



Proof. It is sufficient to prove that H satisfies ()24|) and 



2Re 



Jn 



e3 x — — ) • curio H dx 
0x3) 



8H 



dx? 



dx + 2Re 



Jo, 



2Re 



9H JTa 
-77—H 3 dx 

0x3 

H 3 - — ds. 

0x3 



(27) 



Recall that, for h > sup{x3 : (xi,X2,X3) T G supp(e r — 1)}, there exists a constant < 
i) « 1 such that e r = 1 in (0,2-7r) 2 x (h — rj,h). Hence, a solution H G H^ T (Q,) 3 to the 
problem (fl6|) belongs to ii~p T (fl) 3 n -£fp((0, 27r) 2 x (h — n, h)) 3 due to interior elliptic regularity 
theory. Then one can extend H to a function defined in all of M 3 that is 27r-biperiodic and 
belongs to H^((0, 2vr) 2 x (-00, h)) 3 n # 2 ((0, 2vr) 2 x (h - n, oo)) 3 (This can be seen using [13] 
combined with suitable cut-off arguments.) By abuse of notation, we still denote the extended 
function by H. Let <p G C°°(IR 3 ) be a smooth and non-negative function supported in the 
unit ball and f^ 3 4>dx = 1. For 5 > and x G M 3 let 4> s (x) = 5~ 3 (p(x/5). The convolution 
H s := (f> s * H belongs to H^{Vl) 3 and thus satisfies (|24|) . Then, from Lemma 13.31 and the fact 



that H s -> i? in fl£ T (ft) 3 n # 2 ((0, 2vr) 2 x(h-r), h)) 3 we obtain that 



curl Q (e r 1 curio H s ) curio (e r 1 curl Q H) in L 2 (ri)" 
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Moreover, the convergence in H 2 ((0,2tt) 2 x (h — T],h)) 3 implies that curl Q H s — > curl a H in 
L 2 (r h ) 3 as 6 ->■ 0. Consequently, satisfies (f24"|) . 

It remains to show that -ff also satisfies (|27p . The function if 5 satisfies (|25p and we consider 
the limit of this identity as 5 — > 0. It is easily seen that div a H s div a H = in L 2 (S1). Thus, 
we have 

e 3 -H sS ^° e 3 -H = in ^/ 2 (r ), div a , T fl£ ^° dw a , T H T in i7 p " 1/2 (r ), 

due to the convergence of -ff 5 to H in Hp(Q) 3 . Further, the convergence of H s to H in 
-£fp((0, 27r) 2 x (h — r], h)) 3 and the fact divo-if = on imply that 

— di\ a H° -)■ div a H = - — m H L/z {T h ). 

0x3 0x3 0x3 v 

Plugging in these limits into (|2"5]) shows that ([2T]) holds. □ 



5 Rellich Identity and Solution Estimate 



For establishing uniqueness of solution to the variational problem (|16|) , we derive in this section 
the so-called Rellich identity relating | curl Q H\ 2 and \8H/dx3\ 2 where H is a solution to the 
homogeneous variational problem corresponding to (|16|) . Then, under suitable non-trapping 
and smoothness conditions on the material parameter, integral inequality resulting from this 
identity allow us to obtain estimate for a solution to the homogeneous problem. As mentioned 
in the introduction, the Rellich identity and solution estimate obtained in this section are much 
simpler than the ones in \\.2\ Section 4]. It turns out also that the non-trapping assumptions 
on the parameter material are weaker than the ones in the latter paper. 

The proof of the Rellich identity is based on an integration-by-parts technique that goes 
back to Rellich [15]. Typically, this technique requires more regularity of a solution than just 
to belong to the energy space. In our case we will roughly speaking multiply the Maxwell 
equations (fT7|) . for G = in the right hand side, by X3&H/dx3 and integrate by parts. 

Lemma 5.1 (Rellich Identity). Assume that e^ 1 G Wp'°°(Q) is positive and bounded away 
from zero. Then the following identity holds for all solutions H G if* x(^) 3 t° ^ e homogeneous 
problem corresponding to (|16p . 





2e^ 


dH 


/ 


8x3 


In 





del 



X3- 



dxi 



curl Q if| 2 + 2Re Ve 



dH— 

a — -"3 

0x3 



dx 



+ Re / e 3 x (R a x H) ■ H ds - 2Re / T a (H 3 )H 3 ds = 0. 



(28) 



Proof. Let H G -ffp T (f2) 3 be a solution to the homogeneous problem corresponding to (|16|) . 
First, using integration by parts we have 



f dH T f 

Re / — — • (e3 x curl a H) ds = I 
Jr h dx 3 J Fh 



dH T 



dx? 



ds + Re 



dH T 



r h 9x3 



V Q ,T#3 ds . 
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Note that H satisfies the assumptions of Lemma 14.31 Together with the latter equation we 
obtain 

r Q}J f Q(x £ ~ M f 

2Re / X3— — • curia (e^ 1 curie if) dx = — / — — — — I curl a H\ 2 dx + h / I curl Q H\ 2 ds 
Jn 9x 3 J n 9x 3 J r 



+ 2 / e~ 



dH 



dx? 



9x 3 

dx + 2Re / Ver 1 • — #3 dx - 2Re 
5x 3 



dH 3 

9x 3 



H.ds 



2h 





dH T 


L 


dxs 



ds + 2/iRe 



dH T 

(?X3 



• V a ,T#3 ds . 



We exploit that H solves CCD for G = 0, 

2Re / X3— — • curl a (e7 curl a H) dx = k 2 2Ke / X3 — — ■ H dx = k 2 X3— dx 

Jn 9x 3 J n dx 3 J n 9x 3 



k 2 / \H\ 2 dx +k 2 h / \H\ 2 ds. 



From the last two equations we conclude that 

'd{x 3 e~ l ] 



9X3 



curl Q H\ 2 - k 2 \H\ 2 dx + 2 / e~ l 



dH 



9X3 



a tj 

dx + 2Re / VeJ 1 • — #3 dx 
9x 3 



H 3 ^-ds -2h 

9X3 





dH T 


L 


9x 3 



ds + 2/iRc 



8H 7 

9X3 



• V Q ,r^3 ds 



h [ (I curl a H\ 2 - k 2 \H\ 2 ) ds = 0. 



Due to the variational formulation (|16p for G = 



(e^ 1 ! curl a if| z -k z \H\^)dx + Re / e 3 x (R a x £T) • H ds = (29) 

— 1/2 

since div a ii = in $7 and R a ■ H = in ff p (T/j) due to Lemma 13.31 Adding the last two 
equations yields that the term J n k 2 \H\ 2 dx cancels, and further exploiting dH 3 /dx 3 = T a H 3 
on Th to yields that 



/ 



de r 



X3 ^ — 

n 9x 3 



curl a H\ dx +2 e r 



9F 



9x a 



dx + 2Re / Ver 1 • — ~H~ 3 dx 

9X3 



2Re / T a (H 3 )H 3 ds + Re / e 3 x (R a x H) ■ H ds + 2/iRe 



dH T 

h dx 3 



■ V a , t# 3 ds 



+ (|curl Q F| 2 -A: 2 | J ff| 2 -2 ^ ) 



ds = 0. 
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Recall equality 

| curio, H\ 2 = | curler H | 2 + \cuvl at THs\ 2 + 
Combining the last two equations yields 



dH T 



dx? 



2Rc 



L(H) = h [ 
Jv h 



+ k 2 \H\ 2 - I curU t ^T| 2 - \curl aT H 3 \ 2 )ds 



where L(H) is the left hand side of (|28|) . It remains now to prove that the right hand side 
of the latter equation vanishes. First, we recall from Lemma 13.31 that dH/dx 3 = T a H in 

— 1/2 

H v (Th) which yields that 







L 


dx 3 



I Pn -Hr.nl 

nGA 





oh 3 


L 


dx 3 



^ \Pn.H 3 ^ n \ 
nGA 



Using the latter formulas and replacing k 2 by \n + a\ 2 + /3 2 in the first boundary term in 
yields 

2 



dH T 



dx 3 



+ k 2 \H\ 2 - | curler ff| 2 - Icurla-Ttfgl 2 I ds 



^ [\f3 n H Ttn \ 2 + (|n + a| 2 + /3 2 )(|# T ,n| 2 + l#3,n| 2 ) - |(n + a) x H T>n \ 2 -\n + a\ 2 \H 3>n \ 



nGA 



[(Pn + l/3n| 2 )|^T,n| 2 + |n + a| 2 |tf T ,n| 2 ~ \(n + a) X H T , n \ 2 + /3 2 |tf 3 ,n| 



nGA 



(30) 



On the other hand, due to the divergence-free condition, we have 

[|n + a| 2 |-H Tjn | 2 - \{n + a) x H T) n\ 2 ] = Yl l( n i + «i)#i,n + («2 + a 2 )#2,n| 



nGA 



nGA 



|div QiT i? T || 2 2(rh) = \\dH 3 /dx 3 \\ 2 L2{rh) = ^ |/3„#3,n| 



nGA 



Now substituting the latter equation into (j30|) leads to 



It +k 2 \H\ 2 -\curl atT H\ 2 -\^\ a)T H 3 \ 2 ^ds =2 ^ Pl\H n \\ (31) 

where we exploited that (3 n is either a non- negative or a purely imaginary number. The proof 
is hence finished if we show that X}/3 n >o Pn\H-n\ 2 = (since then L(H) = 0, which is the claim 
of the theorem). First, we compute that 



(e 3 X (R a x H), H) Th = + a ) " H T, n H 3 , n - Y i/3n\H T , n \ 

nGA nGA 

= J2 ^l^3,n| 2 " ^ ^nl^r,n| 2 . 
nGA nGA 
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Since Re (f3 n ) > this implies that 

Im(e 3 x (R a x H),H) Tfl = - ^ Re (M\H n \ 2 < 0, and 

rceA 

Re(e 3 x (R a x H),H) Th = ^ Im (^)|tf 3i „| 2 + ^ Im (/^ |# Tj , 



(32) 
(33) 



neA 



neA 



(The second equation will be exploited later on.) Taking the imaginary part of the variational 
formulation (|16p with G = and F = H, and exploiting Lemma 13.31 we obtain that 



= lm(e 3 x (R a xH),H) r , 



l32l) 



■^Re(/3 n )|tf n | 

neA 



This implies that \H n \ 2 = for all n such that Re (/3 n ) > 0. Since f5 n is either purely imaginary 
or non-negative, we conclude that ^/3 n >o Pn\Hn\ 2 = 0. □ 

The next Poincare-like result is classical (see, e.g., [8j for a proof). 

Lemma 5.2. For «e{«£ Hp(Q) : v\r = 0} there holds 2|| , u||^2(q) < ^ 2 1| du/dx 3 \ \ 2 L 2^y 

The following assumptions on e" 1 will guarantee a stability estimate and a uniqueness 
statement for a solution to the variational problem (|16|) : 



(a) er 1 £ satisfies — < in ft, 



(6) It holds that 



dx 3 

< in a non-empty open ball B C ft, 



(34) 



(c) There exists 5 > 1/2 such that -||Vre7 1 |lL oo (n) 2 ^ — 



de~ l 



dx? 



2 



Remark 5.3. Note that (|34|) f a) implies that e~ l > 1, since, by construction, e~ l = 1 in 
{h — r] < x 3 < h} for some small ij > 0. For the case of periodic non-absorbing structures, 
the main difference between these non-trapping conditions and the ones for the scalar case 
in is the additional condition (|34p (c). This condition arises from estimating the term 
2Re J^Ve" 1 ■dH/dx 3 Jh)dx in the Rellich identity ([28]) using the Poincare-like result above. 
This is natural since the Rellich identity resulting from a similar technique for the scalar 
case [7|/ does not have a corresponding term. 

Let us construct a function e" 1 that satisfies the above assumptions (|34|) . Choose constants 
< h\ < h,2 < h, A > 0, and a C 1 -smooth cut-off function x £ C 1 ((0,27r) 2 ) with compact 
support in (0, 2tt) 2 such that < x < 1 and x = I in (vr/2, 37r/2) 2 . For x = (xi, X2, x 3 ) T £ ft, 
we define 

Xx(xi,x 2 ) + l, 0<x 3 <hi, 
: T 1 (xi,x 2 ,x 3 ) = { A (ffE^f) X(xi,x 2 ) + 1, hi < x 3 < h 2 , 
1, h 2 < x 3 < h. 
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Then e" 1 is a decreasing function that satisfies (|34p (a), and condition (|34|) ( c) is satisfied when 
X > is small enough. Moreover, e~ l also satisfies condition (|34p (b) in (ft/2, 3tt/2) 2 x (hi, fo). 
However, e" 1 does not satisfy the corresponding conditions (7.2)(b,c) in 112$ . which require, 
roughly speaking, strict positivity of de r /dx 3 in (0, 2ir) 2 x (hx,h2) (an arbitrary ball B C O as 
in ()34|) (b) is not sufficient for the proof in \12$ ). 

Lemma 5.4. Assume that e~ l satisfies the three assumptions in f|34j) . Then there exists C > 
(independent of k > 0) such that 

2 



C 



dH 



Ox? 



dx < 



de r 



X3 



dx?. 



curl a H\ dx 



for all solutions H £ i?p T (f2) 3 to the homogeneous problem corresponding to (|16p . 
Proof. We first estimate the two boundary terms in (|28p . We find that 

-2Re / T a (H 3 )lhds = 2 ^ Im (/? n )|ff 3 , n | 2 > 0. 

neA 

Together with ()33p we obtain 

Re(e 3 x (R a x H),H) Th - 2Re / T a (H 3 )~H 3 ds = ^ Im (/? n ,)|# n ,| 2 > 0. 

nGA 



Therefore, from the Rellich identity (|28|) we deduce V(^) < where V(H) is the volumetric 
terms in (|28p . We need now to bound V(H) from below, 



V(H) 



dH 



dx 3 



de r 



X3 



dx ? . 



curl Q #| 2 + 2Re V T £ r ~ 



> 



dH 



dx 3 



del 



X3- 



dx 3 



curl a H\ 2 dx 



dx — 7 



dx 3 



.1 dH T 

dx 3 
2 

- 7 

L 2 (H) 



H 3 + 



de~ l 

r 



de" 1 cff 3 
<9x 3 <9x 3 
2 



dx 



<9x? 



Iff 



L°°(Q) 



^||VT£ r ll|cx»(o)2||-H"3|| 2 :2( n ) - <5 ] 



dH T 



dx 3 



3\\L 2 (n) 



L 2 (fl) 2 



for arbitrary 5, 7 > 0. Poincare's inequality from Lemma 15.21 and the binomial formula imply 
that 



V(H)> 



( 2 ^-|l v T£ r 1 Hi«»(n)a) 



dH, 



dx 3 



+ 



25-1 



8 

dH 3 



dHi 



dx 3 



dx 3 



L 2 {Q.) 



d £ r 1 1 , tti2 1 

x 3 — — I curi a /i I dx 



<9x 3 



dx 



7 



dx 3 



\H 



L°°(Q) 



3Hl 2 (Q) 



Again, we exploit Poincare's inequality, to find that 



7 



dH 3 


2 


de^ 1 


2 


dx 3 


+ 7 

L 2 (f2) 


<9x 3 


L°°(H) 



|ff 3 || L 2 (n) < I 7 + 7y 



5x.s 



dH* 



dx? 
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The minimum of 7 1— > 7 1 + C7 is 2y/C. In consequence, 



V(H)> 



||VT£ r 1 ||ioo( C )2 







/ 




L°°(C) 


Jf7 



<9#3 



dx 3 



+ 



25-1 



8 



8Hq 



dx? 



dx 



dx 



£3 — — I curl a H I dx . 



Finally, assumption (|34[)(c) implies that there exists <5 > 1/2 such that the first bracket on the 
right-hand side is positive. □ 



6 Uniqueness of Solution for All Wave Numbers 

In this section, we prove our main uniqueness result for the electromagnetic scattering prob- 
lem (|16p . under the assumption that e r satisfies (|34p . As mentioned above, corresponding 
uniqueness results that hold for all wave numbers currently exist, to the best of our know- 
ledge, only for absorbing materials, see [16], or simpler two-dimensional structures, see [7]. 

Theorem 6.1. Assume that e" 1 satisfies the assumptions (|34|) . Then problem (|16|) is uniquely 
solvable for all right-hand sides G £ H^(Q,) and for all wave numbers k > 0. 

Proof. Consider a solution H £ ffp T (f]) 3 to the homogeneous problem corresponding to (|16|) . 
Due to Lemma f5.4l and the assumptions on e" 1 we obtain that dH/dx^ = in £1 and curl a H = 
in the ball B (see assumption (|34p (b)). Equation (|17p implies that H vanishes in B, too. 

Since H is independent of x$, it is sufficient to show that H vanishes on T^-^ = 
{(xi,X2, X3) £ : X3 = /i — 77} for some (small) 77 > to conclude that H vanishes entirely in 
Q. If r\ is small enough, then all three components Hj, j = 1, 2, 3, satisfy 

A a Hj + k 2 Hj = 0, A a Hj := Ailj + 2ia • VHj - \a\ 2 Hj, 

in some neighborhood of T^-^. Let us denote by A2 = d 2 /dx 2 + d 2 /dx 2 the two-dimensional 
Laplacian. Since d 2 Hj/dx 2 vanishes, Hj\ r £ H^iT^^) is a weak solution to the two- 
dimensional equation 

A 2 Hj + 2ia ■ V T Hj + (k 2 - \a\ 2 )Hj = on F h _ v , j = 1, 2, 3. 

Standard elliptic regularity results imply that Hj\j, belongs to H^T^^). Moreover, since 
H vanishes in the open ball B and since H is independent of x%, Hj vanish in a non-empty 
relatively open subset of Th- V - 

In this situation, the unique continuation principle stated in Theorem 16.21 (see, e.g., [14] ) 
implies that Hj vanishes on Th-r) for j = 1, 2, 3, and hence H vanishes in Q. □ 

Theorem 6.2. Let O be an open and simply connected set in M 2 , and let u\, ...,u m £ H 2 (0) 
be real-valued such that 

m 

\Auj\ < C^2(\ui\ + \Vu t \) in O for j = l,...,m. (35) 
1=1 



20 



If Uj vanishes in some open and non-empty subset of O for all j = l,...,m, then Uj vanish 
identically in O for all j = 1, m. 
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